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Abstract 



1-^ ■ This paper concerns the global propagation of impact-induced tensile waves in 

"^ ' a kind of phase-transforming materials. It is well-known that the governing system 

"S ; 

jrt I of partial differential equations is hyperbolic-elliptic and the initial-boundary value 

problem is not well-posed at all levels of loading. By making use of fully nonlinear 
stress-strain curve to model this material, Dai and Kong succeeded in constructing 

^ ' a physical solution of the above initial-boundary value problem. For the impact of 

fvj ' intermediate range, they assumed that j3 < 3a in the stress-response function for 



simplicity. In this paper, we revisit the impact problem and consider the propagation 



r~^ . of impact-induced tensile waves for all values of the parameters a and l3. The physical 

^— V ■ solutions for all levels of loading are obtained completely. 
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1 Introduction 

For a thin bar or rod made of a kind of phase-transforming materials, and placed in an 
equilibrium state of uniaxial tension, the curve of the nominal stress and longitudinal strain 
is concave from small to moderate strains but strongly convex for large strains. This kind 
of phase-transforming materials have many applications and have been studied by many 
people (e.g., [I]-[Z] and references therein). It is well-known that the system of conservation 
laws governing the dynamics of the elastic bar is hyperbolic for a standard material for 
which the nominal stress is a monotonically increasing function of the longitudinal strain, 
and is hyperbolic-elliptic for a typical phase transforming material. Abeyaratne et al. [1] 
considered the impact-induced phase transition problem in a semi-infinite bar with a given 
velocity —V at the end. By using tri-linear approximation for the stress-strain curve, they 
constructed a two-wave solution involving a shock wave and a phase boundary. Knowles 
[H] studied the impact-induced tensile waves in a one-dimensional semi-infinite bar made 
of a rubberlike material. Knowles showed that there are three regimes of response related 
to the loading intensity and for the impact of intermediate strength, under the hypotheses 
such as maximally dissipative kinetics or dissipation free kinetics, he selected the unique 
admissible solution among one-parameter family of solutions. 

Recently, Dai and Kong [6] used a fully nonlinear stress-strain curve instead of a 
linearized tri-linear curve in literature, and with the aid of the uniqueness condition on 
phase boundary derived by Dai [5], they succeeded in constructing the unique physical 
solution of the initial- boundary value problem for the loading at all levels. In particular, 
for the intermediate impact velocity, Dai and Kong [5] assumed that the two parameters 
in the stress-response function (see (j2.3p below) satisfy /3 < 3a, and based on this, they 
established the unique physical solution containing a centered rarefaction wave followed 
by a phase-boundary. They also pointed out that the characteristic speeds in the front- 
and back-states of the phase boundary are equal. 

In this paper, we revisit the initial-boundary value problem mentioned above for a 
kind of phase transforming materials. The inequality /? < 3a will be refined and replaced 
by another one and then all physically admissible solutions for all levels of loading are 
obtained completely. Especially, for the intermediate impact, we also construct a physical 
solution involving a centered rarefaction wave followed by a phase boundary. Comparably, 
the geometric character of this phase boundary exhibits a different manner. 



The paper is organized as follows. In Section 2, we recall the initial-boundary value 
problem modeling the impact-induced tensile waves in a semi-infinite bar made of a kind 
of phase-transforming material. Section 3 is devoted to the establishment of admissible 
solutions for all loading at all levels. Particularly, in the intermediate case, we completely 
construct physical solutions for all parameters in the stress-response function. We conclude 
the paper in section 4 and some related remarks are also given. 

2 The initial-boundary value problem 

Consider the time-dependent deformation of a one-dimensional elastic bar from the natural 
state which occupies the nonnegative j;-axis. At the time t, a particle at x in the reference 
state is carried out to a; + u{t,x), where u denotes the longitudinal displacement of the 
particle. Hence the strain 7(t, x) and the particle velocity v{t, x) are defined by 7 = Ux 
and V = ut respectively. To ensure the deformation is one-to-one, we suppose 7 > — 1. 
Let a{t, x) stands for the nominal stress at the time t for this particle. In Lagrangian 
description, the governing system of equations of motion is given by 

PVt = CTx, It = Vx, (2.1) 

where 7, v are smooth, and the jump conditions 

s[7] + [v] = 0, sp[v] + [a] = (2.2) 

at a moving strain discontinuity whose referential location is x = s{t) at time i., in which 
[/] = f{t,s{t) -|- 0) — f{t,s{t) — 0), p is the constant density and s = s'{t) is denoted by 
the Lagrangian velocity of the discontinuity. 

The material of the bar is taken to be elastic so that a = (7(7), where the stress- 
response relation a = (7(7) is assumed to be a C^-smooth function. Thus, the system 
(|2.ip is hyperbolic-elliptic for typical phase-transforming materials for which the stress- 
strain curve generally has a peak-valley combination. 

As in Dai and Kong [6], this paper concerns the following nonlinear stress-strain rela- 
tion in order to model two phase materials 

a{^) = E(^-l{a + l3)j^ + al3j], (2.3) 
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where E is the Young's modulus for infinitesimal strains, a and f3 are constants satisfying 

< a < /? < 00. (2.4) 



a = (7(7) 



Fig. 1: The stress-strain curve (7(7) 

In general, some parts of the curve may lie below the strain axis. The restrictions of the 
stress-response function ()2.3p to the intervals [0, a] and [/3, 00) are called a-branch and 
/3-branch respectively. Obviously, 17(7) is concave for — 1 < 7 < (a + /3)/2 and convex for 
7 > (a + /3)/2. see Fig.l. 

The sound-wave speed 0(7) for the bar characterized by ()2.3p is given by 

c(7) = ^/a'{-,)/p = coVh - «)(7 - /3)/a/3 for 7 G [0, a] U [/3, 00), (2.5) 

where cq is the speed of small amplitude waves at the undeformed state 



Co = c 



(0) = yjEaP/p. 



As in [6], for the system (j2.ip . the impact problem corresponds to the following initial- 
boundary conditions 

7(0,3;) =0, i;(0,3;) =0 forx>0, (2.6) 



t,(t,0) = -y fort>0. 



(2.7) 



where V is the impact velocity and assumed throughout to be positive. We shall seek the 
solution of the initial boundary value problem (|2.ip . (|2.6p - (|2.7p on the first quadrant. 



3 The impact-induced tensile waves 

In this section, we exhibit the structure of solutions of the impact problem ()2.ip . ()2.6p - ()2.7p 
for all values of a, /3 in the stress-response function (|2.3p under the assumption (|2.4p . 



Introduce 



K = /" c(7)d7 = ci f VC^ - WC, K* = 2cihJ^ - 1, K„ = Scih^h^ - 1, (3.1) 



where 



A (/3-a)^ ^ /3 + a 

ci = — ^=-co, /i = - — -. (3.2) 

4:^/ap p — a 



Remark 3.1 The values of constants V^,, y*^, and Ki,*^, coincide with those in f^. 

It is easy to derive the fohowing lemma. 

Lemma 3.1 Under the assumption ^2.4\ ), Ki.,Kk** are positive constants and satisfy 

< K < 14„ < oo; (3.3) 

if furthermore we assume that 

/5<^.. (3.4) 

where /i* G (\/3, 2) is the unique solution to V^ = Kk*, then T4* makes sense and it holds 
that 

< K < K* < K** < oo. (3.5) 

Although for weak impact (^ ^ V*) and strong impact V ^ V***, the solutions for the 
impact problem are the same as constructed in Dai and Kong [6|, we still state them in 
the sequel for completeness. 

3.1 Weak impacts 

As in [6], the solution contains a centered rarefaction wave and is given by 



{v{t,x),'y{t,x)) = < 



(-y,7i) forO^C^ei, 

[(0,0) fore ^6, 

where ^ = x/t, ^2 = cq and v, 7, 71, ^i are defined by 

/■CO rii 

cmO) = e, m = / Cl'iOdC, / c(7)d7 = V and 71 = 7(6). (3.7) 

J( Jo 

The interval (0, K] is the regime of weak impacts. As the impact velocity V £ (0, K], the 

strain 71 satisfies 71 ^ a. 



3.2 Strong impact 

As in [6], the solution contains a shock wave and is given by 

(-y,7~) for0^x<st, 



{v{t,x),-f{t,x)) 



(0,0) 



(3.8) 



for X > st, 



where 7 and s are defined by 

7"co-\/7^V3 - (a + /3)7~/2 + 0^5 = y^^ 



and 



7" 



(3.9) 



Obviously, this pure-shock- wave solutions are constructed for which 7^ > |(q; + /?) and 
therefore V > Ki,**. Moreover, the shock wave is supersonic in the Lagrangian sense with 
respect to the undisturbed state ahead of it, i.e. cq < s < 0(7"). In particular, when 
7~ = |(a + /?), the impact velocity V = K^,* and the discontinuity x = st is called a 
degenerate shock wave since s = cq and s < 0(7^ ). 



3.3 Impact of intermediate strength 

In this subsection, we turn to construct the solution of the initial-boundary value problem 
([2rT]) . (fT6]) - (f2Tl) to fill the gap of impact velocity V^ < V < K**- We now attempt to 
construct one solution involving a centered rarefaction wave followed by a phase-boundary, 
as shown in Fig. 2. Precisely speaking, we propose the following form of solution 



{v{t,x),-f{t,x)) 



for ^ X < st, 
for st < X ^ ^li, 
for ^it ^ x ^ ^2*, 
for X ^ ^2^) 



(«+,7+) 

,(0,0) 
where 7", s, ^1, u and 7 are to be determined and, by continuity, ^2 = co,7~' 

Assuming that 7 and v are continuous everywhere except across the phase boundary, 
we have ^1 = 0(7^^). As in the discussion of the centered rarefaction wave arising from 
weak impact [^, we obtain that 



(3.10) 



7(6) and 



v{^i 



CO 



C7'(C)^C 



c(7)d7. 



(3.11) 
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Since the phase boundary is to trail the centered rarefaction wave, we require that 



< s ^ 6 = c(7+). 



(3.12) 



phase boundary 
^x = st 



X = ^it 

rarefaction wave 




a- = 6i 



Fig. 2: Solution with a centered rarefaction wave followed by a phase boundary 



Using ()3.1ip in the jump conditions ()2.2p applied at the phase boundary x = st leads to 



(7+ -7 )s 



^(7)^7 — y 



0. 



(3.13) 



a(7+) - a(7-) - ps { I c{^)dj - ^ 5> = 



(3.14) 



These are two equations for the front- and back-state strains 7 at the phase boundary; 
they involve the phase boundary speed s as an unknown function of V, while the impact 
velocity y as a given datum. 

Eliminating the common contents of the braces between ()3.13p and ()3.14p yields the 
standard formula relating the phase boundary speed and the slope of the chord connecting 
the two points on the stress-strain curve that correspond to the front- and back-states of 
the phase boundary 

ps' = <^^l-<^~\ (3.15) 

7+ -7 

It is easy to see that the equations (|3.13p - (|3.14p are equivalent to either (|3.13p or (|3.14p 
together with ()3.15p . 

For the material governed by (|2.3p . p.l3p and (|3.15p take the respective forms 

(7+-7-)s+/ c(7)d7 + y - 14 = 0, (3.16) 



+ 7 



a + 13 



3a(3 



s^ + 3 



/? — a 



(3.17) 



We shall view (l3l^ and ([317)1 as the system (f3T^ and (f3H]l . 



For each V between y* and 14**, we shall show that there is a one-parameter family 
of solutions to the initial-boundary value problem ()2.ip . ()2.6p - ()2.7p . So a mechanism is 
needed to select one physical unique solution. Abeyaratne and Knowles introduced the 
concept of driving force, which is defined via the dissipation rate (cf. [2]-Il]); and the 
kinetic relations to select the unique admissible solution. The dissipation rate can be 
written as 

D{t) = f{t)s'{t), (3.18) 

where /(t) is called the driving force per unit cross-sectional area acting at time t on the 
moving strain discontinuity. For a phase-transforming material governed by ()2.3p . direct 
computation shows that / is given by 

/ = §(7+ - 7")' {(« + /3) - (7+ +7")} • (3.19) 

Since s > 0, the physical admissibility in the sense of Knowles (see for definition) 
requires that D{t) ^ 0. Since in (|3.19p one has 7" > 7"^, it implies that 

7++7-^a + /3. (3.20) 

If in (|3.12p one represents s through (|3.17p and c{j'^) through (j2.5p . one finds, after 
simplification, the further restriction 

27++7-^^(a + /3). (3.21) 

Moreover, in order to keep s > in (|3.17p . direct computation gives the following restric- 
tion 



7 > air) = -^— + 2 



2i±^_,+ +/3(M^-,+V/!^ + ,H 



(3.22) 



Pairs (7~'^,7~) satisfying the inequalities (|3.20p - (j3.22p as well as 7"*^ ^ a,7~ ^ /? corre- 
spond to points in the quadrilateral BCEF, the triangles BCE and BCF in the 7^, 7~-plane 
for three cases respectively, see Fig 3. 

The quadrilateral and triangular regions in Fig. 3, along with the portions of straight 
lines that represent the centered rarefaction wave and pure shock wave solutions, can be 
mapped through (|3.16p and (j3.17p into a plane in which the natural Cartesian coordinates 
are the dimensionless phase boundary velocity s/c2 and the dimensionless impact velocity 




O 7* " 



(a) /3 < (2 + V3)a 



7+ O tt 



(b) I3={2 + V3)a 



7+ O a 



(c) /3 > (2 + v^)a 



Fig. 3: Pairs (7''',7~) satisfying i3.13\) - (3.14\ l must respectively correspond to points in the 
quadrilateral BCEF, the triangles BCE and BCF in the ^^^^^ -plane for the cases (a), (h) and 
(c). Points of the form (0,7^) with 7~ ^ |(q; + (3), lying on the vertical axis above the point 
C, correspond to pure shock wave solution. Points on the line OA: 7^ = 7" with < 7'*' < a 
correspond to the centered rarefaction wave. The curve BE or BE is given by 7^ = .9(7^)- 



V/ci, where C2 = 2ci/(/3 — a). Consider first the centered rarefaction wave solutions. 
Setting 7+ = 7~ = ^^j£ — ^^i] in ()3.16p and ()3.17p gives rise to the parametric equations 



C2 



V^F^, 



I Cl Jr, 



1 s; ?7 ^ /i. 



(3.23) 



Next we map the segment BC for which 7 = |(a + /3) — 27"*" with < 7^ ^ a. By 



(J3.16p - (j3.17p . the image of this segment in the dimensionless s,V-plane is the following 
curve 



C2 

V_ 

(. Cl 



V^f^, 



3?7^r/2 - 1 + / ^C^ - 1 dC, 
J-n 



li^i]i^h. 



(3.24) 



From p.l6p - p.l7p . the image of segment CD for which 7"*^ = 0,7 ^ |(a + /3) is found 
to be the curve represented by 



C2 V 3 



— = {ip + h)-, 
Cl C2 



(f ^ 2/1, 



(3.25) 



where 



"P 



27- - (g + /3) 
(3 — a 



(3.26) 



By the same method, the segments CE in the Fig 3. (a) and (b) can be mapped into the 



following parameterized curve 



C2 

- = {ip + h)-, 
Cl C2 



(^2 _ /j^ + /j2) _ 1^ 



h^ip s^2h. 



(3.27) 



Form (j3.16p . (j3.17p and (|3.22p . the segment CF in Fig 3. (c) may now be mapped into 
the curve in the dimensionless s, T^-plane 



— = {ip + h)-, 
Cl C2 



[h + v^3(4 - /i2)] ^ 99 ^ 2/i. 



(3.28) 



In Fig 3. (a), the image of segment EF now takes the following form 



s 

C2 

v_ 

k Cl 



1, 



2r/ 



r?- 



1+ / VC^-i^, 

Jv 



(3.29) 



Finally, the curve BF or BE in Fig 3. can be represented in the following parametric form 

(3.30) 



f s 

C2 



K Cl L 



The curves BF and BE in Fig 3. (a)-(b), BF in Fig 3. (c) are parameterized through 
p.30p with T] G [1, -v/3] and [l,/i] respectively. 

The following Fig. 4 and Fig. 5 show the image just constructed in the dimensionless 
s,y-plane of the locus in the 7"*",7~-plane describing the various types of solutions as 
presented in Fig 3. 

It is worthy pointing out that if, for given s and V, there is a physically admissible 
solution 7''",7~ of the equations (|3.16p . (j3.17p with ^ 7^ ^ a and for which < 
s ^ c(7+), then it has been shown that the point corresponding to s,V must lie in the 
quadrilaterals B'C'E'F' or the triangular regions B'C'E', B'C'F' in Fig. 4 and Fig. 5. 
While, the converse is not proved. This means that for every pair s, V corresponding to a 
point in the above regions, there is an acceptable solution 7^,7" of ()3.16p . (j3.17p . 

Now we refine the inequality /3 < 3a of Dai and Kong [6j into (|3.4p . In other words, 
under the assumption /3 < -j^ji-j-a, by means of the same method, we can also construct 
the unique physical solution involving a centered rarefaction wave followed by a phase 
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Fig. 4: Pairs (s/c2,V^/ci) corresponding to points (7^,7 ) m the locus of Fig. 3 (a). Cases (ai), 
(a2) and (a^) correspond to V^ < T4*, V^ — 14* and V^ > 14* respectively. 
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(b') /3 = (2 + ^)a 



(c') I3>(2 + V3)q 



Fig. 5: Pairs {s/ C2,V/ci) corresponding to points (7''', 7 ) in the locus of Fig. 3 (b) and (c). In 
the present cases, the impact velocity T4* does not make sense. 

boundary, which moreover satisfies the following unique condition on the phase boundary 
derived by Dai [5] 

7++7- = a + /3. (3.31) 

This solution corresponds to the arc GE' in Fig 4. (ai), while the curve O'A' corresponds 
to the centered rarefaction wave for weak impact; the curves C'E' and C'D' relate to the 
admissible solution and pure shock wave solution presented in [6]. Here we would like to 
point out that the mixed centered rarefaction wave and phase boundary solution presented 
in [6] in fact satisfies the dissipation-free kinetics introduced by Abeyaratne and Knowles 
(cf. [2]-[l] or [8]), i.e., the driving force / = 0. 

In the reminder of this subsection, we aim to construct the unique physical solution 
to ([2J]) . ([2:6]) and ([S^D for the left case f3 ^ j^a. For the impact of intermediate 
strength, there is also a one-parameter family of two- wave solutions (see Fig. 4 and Fig 
.5). In order to select the unique admissible solution, we apply the maximally dissipative 
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kinetics, which implies that the phase boundary speed should satisfy (cf. [8]) 

s = 0(7"^). 



(3.32) 



The equation (|3.32p means that the phase boundary in fact coincides with the tail of the 
centered rarefaction wave. By (j3.15p and (j3.32p . we have the following important identity 



27+ + 7" = -(a + /3). 



(3.33) 



Next we shall show that the equations (j3.15p and p.l6p have a solution (7+ , 7 ) 
satisfying 7+ G (0, a). Making use of ()3.32p and ()3.33p . it follows from ()3.16p that 



-(a + /3)-37^ 



c(7+) 



c(7)(i7 + y-T4. 



(3.34) 



Let 



Gil' 



-(a + /3)-37^ 



-il'^) - / c(7)d7 



V.-V, 



it is easy to observe that G{'y^) is a monotonically decreasing function as 7"*" G (0, a). 
Moreover, we have 

G{a) =V^-V <0, G(0) = ^(a + P)co -V = V^^^-V>0. 

Thus, it is found that there is a unique 7"^ G (0, a) such that (|3.34p holds. Once 7"*^ has 

been solved, 7^ can be determined from (j3.33p immediately. Since /i* G (\/3, 2) in Lemma 

3.1, we have 

h. -\-^ 

(3.35) 






K - I 
Hence, noting < 7+ ^ a and (j3.35p . it follows from (|3.33p that 

3/ ns _ 3/3 — a 

-(a + /3)>7 > — - — >a + /3. 



(3.36) 



which implies 

7+ + 7^>a + /3. (3.37) 

Therefore, the dissipation rate D{t) > and this means that for the case /3 ^ ^*]^^ a, the 
constructed unique solution (j3.10p is physically admissible. 

On the other hand, direct calculation shows that 0(7"*") = s < 0(7"), which implies 
that the phase boundary speed is equal to the characteristic speed in the front-state 
and subsonic in the back-state of the phase boundary. This fact can be viewed as the 
geometrical feature of the phase boundary under the assumption /? ^ if^a, while for 
the case /? < jf^a, the phase boundary exhibits a different manner, i.e., 0(7^) > s and 
c(7+) = c(7-)(cf. i). 
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Remark 3.2 Using the method given in |^, we also can investigate the global structure 
stability of the physical solution constructed in this section. The details are omitted for 
space limitations. 

4 Concluding remarks 

This paper concerns the propagation of impact-induced tensile waves in a kind of phase 
transforming materials. For the impact velocity of intermediate strength V^, < V < T4**, 
the biggest problem under consideration is how to select the unique physical solution. 

For the material governed by (|2.3p . we consider the impact problem by dividing the 
parameters in stress-response function into two cases: /3 < -jf^a and /3 ^ h-i '^- ^^ both 
cases, we can establish the mixed centered rarefaction wave and phase boundary solution 
()3.10p . It is pointed out the geometric features of the two phase boundaries are different. 
The mixed-type solution in the former case indeed satisfies the dissipation free kinetics 
and the latter one corresponds to the maximally dissipative kinetics. Precisely speaking, 
direct computation shows that this kinetic relation is given by 

212 



■' 4 



1+' ' 

■ C2 



Most recently, Knowles [9] studies the maximally dissipation kinetics for general nonlinear 
elastic bar and points out it depends strongly on the qualitative nature of the stress-strain 
relation. 

Finally, we would like to mention that for the intermediate impact, under the maxi- 
mally dissipative kinetics hypothesis, the solution containing a centered rarefaction wave 
and a phase boundary can be constructed for all values of a, f3, see the curves O'B'C' in 
Fig. 4 and Fig. 5; while under the dissipation free kinetics, construction of this mixed-type 
solution is not obvious except the case /? < j^^^Oi as discussed in the last section. 
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